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A random matrix model to describe the coupling of m-fold symmetry is con- 
structed. The particular threefold case is used to analyze data on eigenfrequencies 
of elastomechanical vibration of an anisotropic quartz block. It is suggested that 
such experimental/theoretical study may supply powerful means to discern intrinsic 
symmetries of physical systems. 

The standard ensembles of Random Matrix Theory (RMT) [l| have had wide appli- 
cation in the description of the statistical properties of eigenvalues and eigenf unctions of 
complex many-body systems. Other ensembles have also been introduced 0], in order 
to cover situations that depart from universality classes of RMT. One such class of en- 
sembles is the so-called Deformed Gaussian Orthogonal Ensemble (DGOE) 0, 0, 0, 0] 
that proved to be particularly useful when one wants to study the breaking of a discrete 
symmetry in a many-body system such as the atomic nucleus. 

In fact, the use of spectral statistics as a probe of symmetries in physical systems 
has been a subject of intensive experimental and theoretical investigation following the 
pioneering work of Bohigas, Giannoni and Schmit |7| which showed that the quantal 
behaviour of classically chaotic systems exhibits the predictions supplied by the RMT. 
Examples of symmetry breaking in physical systems that have been studied include nuclei 
@,@], atoms and mesoscopic devices such as quantum dots [l2| . 



In the case of nuclei, the Mitchell group at the Triangle Universities Nuclear Labo- 
ratory j^, loj , studied the effect of isospin symmetry breaking, in odd-odd nuclei such as 
"^^Al. They detected the breakdown of this important symmetry by the applications of 
two statistics: the short-range, nearest neighbor level spacing distribution (NND) and 
the long range Dyson's A-statistics jl, [9|. These results were well described by a DGOE 
in which a pair of diagonal blocks is coupled. The strength of the coupling needed to 
account for the symmetry breaking can be traced to the average matrix element of the 
Coulomb interaction responsible for this discrete symmetry breaking 0, 13]. The justi- 



fication for the use of block matrices to describe the statistics of a superposition of R 
spectra with different values of the conserved quatum number can be traced to Refs. 

Ijj. In the case of non-interacting spectra, i.e. if the quantum number is exactly 
conserved, the answer is a superposition of the R spectra. Since the level repulsion is 
present in each one of the R spectra, their superposition does not show this feature. 
Thus, we can say that for each spectra of states of a given value of the quantum number, 
one attaches a random matrix (GOE). For R spectra each of which has a given value of 
the conserved quantum number, one would have an R x R block diagonal matrix. Each 
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block matrix will have a dimension dictated by the number state of that spectra. If the 
quantum number is not conserved then the R x R block matrix acquires non-diagonal 
matrices that measure the degree of the breaking of the associated symmetry. This idea 
was employed by Guhr and Weidenmiiller [13] and Hussein and Pato [1] to discuss isopin 
violation in the nucleus '^^Al. In reference [3|, the random block matrix model was called 
the Deformed Gaussian Orthogonal Ensemble (DGOE). 

In order to study transitions amongst universal classes of ensembles such as order- 
chaos (Poisson— i>GOE), symmetry violation transitions (2G0E— i>lGOE), experiments 
on physical systems are more complicated due to the difficulty of tuning the interaction 
(except, e.g. in highly excited atoms where the application of a magnetic field allows 
the study of GOE-GUE transitions). To simulate the microscopic physical systems, one 
relies on analog computers such as microwave cavities, pioneered by A. Richter and 
collaborators 15| and acoustic resonators of EUegaard and collaborators [l^, 18]. It 



is worth mentioning at this point that the first to draw attention to the applicability of 
RMT to accoustic waves in physical system was Weaver [l9| . 



In the experiment of EUegaard et al. [17| what was measured were eigenfrequencies 



of the elastomechanical vibrations of an anisotropic crystal block with a D3 point-group 
symmetry. The rectangular crystal block employed by Ellegard was so prepared as to 
have only a two-fold flip symmetry retained. Then, to all effects, the quartz specimen 
resembles a system of two three-dimensional Sinai billiards. The statistical treatment 
of the eigenfrequencies of such a block would follow that of the superposition of two 
uncoupled GOE's. 

Then, by removing octants of progressively larger radius from a corner of the crystal 
block this remnant two-fold symmetry was gradually broken. The spectral statistics show 
a transition towards fully a chaotic system as the octant radius increases. What was then 
seen was that the measured NND is compatible with a two block DGOE description but 
the A-statistics was discrepant. This discrepancy was attributed to pseudo integrable 
behavior and this explanation was later implemented with the result that the long-range 
behavior was fitted at the cost, however, of loosing the previous agreement shown by the 
NND (23,]. 

Here we reanalyse this experiment following the simpler idea of extending the DGOE 
matrix model [5| to consider the coupling of three instead of two GOE's . We show that 
within this extension both, the short- and the long-range statistics, are reasonably fitted 
suggesting that the assumption of the reduction of the complex symmetries of anisotropic 
quartz block may not be correct. Our findings have the potential of supplying very precise 
means of testing details of symmetry breaking in pysical systems. 

To define the ensembles of random matrices we are going to work with, we recall 
the construction based on the Maximum Entropy Principle [31], that leads to a random 
Hamiltonian which can be cast into the form 

H = Ho + XHu (1) 

where the block diagonal Hq is a matrix made of m uncoupled GOE blocks and A (0 < 
A < 1) is the parameter that controls the coupling among the blocks represented by the 
Hi off-diagonal blocks. For A = 1, the Hi part completes the Hq part and H = H^^^ . 

These two matrices Hq and Hi are better expressed introducing the following m pro- 
jection operators 

P^ = Y,\J><]l (2) 

3&h 



where /j defines the domain of variation of the row and column indexes associated with 
ith diagonal block of size Mj. Since we are specifically interested in the transition from 
a set of m uncoupled GOE's to a single GOE, we use the above projectors to generalize 
our previous model jl, 0] by writing 



GOE 



Pi 



(3) 
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and 
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where Qj = 1 — Pj. It is easily verified that H = H'^^^ for A = 1. 

The joint probability distribution of matrix elements can be put in the form j^, 20 



(5) 



(6) 



P{H, a, (3) = Z^^ exp {-atrH^ - (3trHl) 
with the parameter A being given in terms of a and (3 by 

X = {l + ^)-^/\ 
a 

Statistical measures of the completely uncoupled m blocks have been derived. They 
show that level repulsion disappears which can be understood since eigenvalues from 
different blocks behave independently. In fact, as m increases the Poisson statistics are 
gradually approached. In the interpolating situation of partial coupling, some approxi- 
mate analytical results have been derived. In Ref. [13], for instance, it has been found 
that the density p{E) for arbitrary A and m is given by 
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is Wigner's semi-circle law with a = \J N/a and 
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Eq. ([5|_^can be used to calculate exactly analytically the NND for 2x2 and 3x3 
matrices [6|. For the 2x2 case the DGOE, Eq. ([5]), gives 
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where /q is the modified Bessel function, whose asymptotic form is 



Thus,there is no level repulsion for — > oo, P2x2('5, oo) = |aexp (— fs^), which repre- 
sents the 2x2 Poisson distribution where the usual exponential is replaced by a Gaussian. 
The prefactor is just 1 if 2a is taken to be tt. In the opposite limit, /3 — > 0, Iq{x) ~ 1 — 
and one obtains the Wigner distribution, 

P2x2(s,/3 0) ^ ^exp {-^s^) s (12) 

Note that the parameter A of eq ([6]) is if /5 is oo and 1 if /5 is 0. 

For higher dimensions Eq. ([5]) can only be used for numerical simulations. This is what 
we are now reporting, using 2 and 3 bolck matrices of sizes 105 x 105 and 70 x 70 each, 
respectively. The size of the whole matrix is 210 x 210. Further, we take an ensemble of 
1000 elements and fix a to be 1. We apply our model to analyse the eigenfrequency data 
of the elastomechanical vibrations of an anisotropic quartz block used in [l7| . 

In this reference in order to break the flip symmetry of the crystal block gradually 
they removed an octant of a sphere of varying size at one of the corners. The rectangular 
quartz block has the dimensions 14 x 25 x 40 mm^. The radii of the spheres containing 
the octants are r = 0.0, 0.5, 0.8, 1.1, 1.4 and 1.1 mm representing figures (a) — (/). Figs. 
Ix and 2x of Ref. [t?] correspond to an octant of a huge sphere of radius r = 10.0 mm, 
whose center is inside the crystal and close to one of the corners. They found 1424, 1414, 
1424, 1414, 1424 and 1419 frequency eigenmodes in cases (a) — (/), respectively. These 
eigenfrequecies were measured in the frequency range between 600 and 900 kHz. Thus 
the average spacing between the modes is about 214Hz. The histograms and circles in 



the two figures of Ref. [iTj represent the short-range nearest-neighbor distributions P{s) 
(Fig. 1) and the long range A3(L) statistics (Fig. 2). In our DGOE simulation the 
unfolding of the calculated spectra is performed with the DGOE density given by Eq. 
(7) above. 

In figures 1 and 2, we show the results of our simlulations as compared to the data 
of Ellegaard et al. [l7|] for the spacing distribution and in figures 3 and 4 the long 
range correlation exemplified by the spectral rigidity A3(L). We simulate the gradual 
breaking of the 2- or 3-fold symmetry by changing the value of the parameter A above. 
We see clearly that in so far as the ^^{L) is concerned a 3-GOE description works much 
better than a 2-GOE one. It is clear, however that both descriptions fall below the data, 
specially at large L. We shall analyse this discrepancy in the following using the missing 
level effect UK. 

It is often the case that there are some missing levels in the statistical sample analysed. 
Such a situation was addressed recently by Bohigas and Pato [2l| who have shown that 



if g fraction of the levels or eigenfrequencies is missing, the A3(L) becomes 

^l{L)=g^^ + {l-gf^,(^-^^. (13) 

The presence of the linear term, even if small, could explain the large L behavior of the 
measured A3(L). We call this effect the Missing Level (ML) effect. Another pos sible 
deviation of A3 could arise from the presence of pseudo-integrable effect (PI) [22|, |23l . 
This also modifies A3 by adding a Poisson term just like Eq. ( |T3i) . In the following we 
show that there is no need for the PI effects to explain the large-L data on the A3 if the 
ML effect is taken into account. 

We take a study case Figs. 3b and 4b which correspond to r = 0.5 mm and where 1414 
frequency eigenvalues were found. We consider this a potential ML case and take for A3, 



the expression given in Eq. fll3l) and apply to our simulations. We find perfect fit to the 
data, if g is taken to be 0.1, namely only 90% of the eigenfrequencies were in fact taken 
into account in the statistical analysis. There is, threfore, room to account much better 
for all cases (Fig. 2a, 2c, . . . ) by appropiately choosing the correponding value of g. We 
have also verified that if a 2G0E description is used, namely, m = 2 , then an account of 
the large-L behaviour of A3 can also be obtained if a much larger number of levels were 
missing in the sample. In our particular case of Fig. 2b, we obtained g = 0.3. This is 3 
times larger than the ML needed in the 3G0E description. We consider the large value 
of g needed in the 2G0E description, much too large to conform to the reported data in 
Ref [13]. Figure 5 summarizes our the above. 

It is therefore clear that the 3G0E description of the spectral rigidity of the eifen- 
frequency spectra of [17] for the crystal block does work very well if a small fraction of 
the levels is taken to be missing, without resort to pseudointegrable trajectories or levels 
that do not feel the symmetry breaking [23|. On the other hand, the 2G0E description, 
which does as good as the 3G0E one in fitting the measured P{s), fails dramatically in 
accounting for the spectral rigidity, even if as much as 30 per cent of the levels are taken 
as missing. 

In conclusion, a random matrix model to describe the coupling of m-fold symmetry is 
constructed. The particular threefold case is used to analyse data on eigenfrequencies of 
elastomechanical vibration of a anisotropic quartz block. By properly taking into account 
the ML effect we have shown that the quartz block could very well be described by 3 
uncoupled GOE's , which are gradually coupled by the breaking of the three-fold sym- 
metry (through the gradual removal of octants of increasing sizes), till a IGOE situation 
is attained. This, therefore, indicates that the unperturbed quartz block may possess 
another symmetry, besides the flip one. A preliminary version of the formal aspect of 
this work has previously appeared in [2J . 
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FIG. 1: Nearest Neighbour Distributions. Histograms show data (a)-(x) from Ref. [l7[. Thick 
histograms show the three coupled GOE fits to the data carried out using the DGOE numerical 
simulations using Eq. ([5]). Also shown as the full thin line the three uncoulped GOE P{s). 
In graph (x) the dotted line is the Poisson distribution, the dashed line is the two uncoupled 
GOE -P(s). The very thin line is Wigner distribution which is hidden behind histograms. The 
values of A that adjust the data are 0.0032, 0.0071, 0.0158, 0.0250, 0.0333, 0.9950, 1.000 for cases 
(a)-(x). See text for details. 




FIG. 2: Nearest Neighbour Distributions. Histograms show data (a)-(x) from Ref. [l7[. Thick 
histograms show the two coupled GOE fits to the data carried out using the DGOE numerical 
simulations using Eq. ([5]). Also shown as the full thin line the two uncoulped GOE P{s). In 
graph (x) the dotted line is the Poisson distribution, the dashed line is three uncoupled GOE 
P{s). The very thin line is Wigner distribution which is hidden behind histograms. The values 
of A that adjust the data are 0.000, 0.0258, 0.0200, 0.0400, 0.0705, 0.0600, 1.000 for cases (a)-(x). 
See text for details. 




FIG. 3: Spectral Rigidities. The thick hnes are the DGOE simulation for the three coupled 
GOE's.The same values of A as in Fig. 1 were used. The thin lines correspond the three 
uncoupled GOE's case. The data points are from Ref. [l7|. See text for details. 




FIG. 4: Spectral Rigidities. The thick hnes are the DGOE simulation for the two coupled 
GOE's. The same values of A as in Fig. 2 were used. The thin lines correspond the two 
uncoupled GOE's case. The data points are from Ref. [l7[. See text for details. 
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FIG. 5: The ML effect. The data points correspond to case (b) of Ref. 17|, r = 0.5 mm. The 
full line corresponds to our three coupled GOE's fit with A = 0.0071, figure 3b and g = 0.10. 
The dashed line corresponds to our two coupled GOE's fit with A = 0.0258, figure 4b and 
g = 0.30. See text for details. 



